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Abstract 

We show that polarons can sustain loop-like configurations in flex- 
ible bio-polymers and that the size of the loops depend on both the 
flexural rigidity of the polymer and the electron-phonon coupling con- 
stant. In particular we show that for single stranded DNA (ssDNA) 
such loops can have as little as 10 base pairs. For polyacetylene the 
shortest loop must have at least 12 nodes. We also show that these 
configurations are very stable under thermal fluctuations and can fa- 
cilitate the formation of hairpin-loops of ssDNA. 

1 Introduction 

Conformational transitions of biopolymers as a result of the coupling be- 
tween the electronic and elastic degrees of freedom are important for un- 
derstanding native states of globular proteins and secondary structures of 
biopolymers such as DNA and RNA. In an attempt to understand toroidal 
states of DNA the globule-coil transition for semi-flexible polymers in poor 
solvents has been explored using Brownian dynamics simulations [1, 2]. The 
intermediate states arising in these systems have also been classified [1, 2]. 
However the collapse transition in polymers induced by polarons has been 
less explored [3]. 

Polarons are the result on the interaction between a free electron in the 
conducting band of a polymer chain and the phonons of that chain. They 
were discovered by Davydov [4, 5, 6] who proposed them as a mechanism to 
explain how energy can be transported along alpha- helices in living cells. 
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In this paper we explore the possibility of polaron induced polymer loop 
formation and stabilisation arising in semi-flexible chains. Our model is sim- 
ilar to the model proposed by Mingaleev et al. [3] who generalised the orig- 
inal model of Davydov [4, 5] by incorporating long ranged electron-phonon 
interactions. In their work Mingaleev et al. showed that at zero temper- 
ature polarons can induce a spontaneous bend in a straight chain if the 
bending modulus is less than a critical threshold. A careful examination of 
the model however reveals that realistic polymers e.g. DNA and polyacety- 
lene are more rigid having their bending modulus twice and twenty times 
above the threshold value respectively. Therefore though interesting from a 
theoretical point of view, the spontaneous bending on polymers induced by 
polarons is limited in scope when applied to physical systems. 

In a recent paper [7] we showed that the Mingaleev et al. model can 
explain spontaneous polaron transport on a chain having a bending gradi- 
ent, e.g. alpha-helices of light harvesting proteins. In this case, the bending 
of the chain is generated by the natural folding of the protein which can 
induce a spontaneous polaron displacement. We showed that with the poly- 
mer configuration frozen in, the polaron spontaneously accelerates along the 
bending gradient, and gets reflected across sharply kinked junctions. Fur- 
ther we showed that at finite temperatures the polaron undergoes a biased 
random walk to a region of high curvature. 

While polarons are not able to induce spontaneous conformational tran- 
sitions in DNA and polyacetylene, on account of their rigidity, they might 
sustain a folded configuration that might have been formed by other means 
e.g. thermal fluctuations, or mechanical stress. This is particularly true for 
ssDNA whose bending modulus is only twice as large as the threshold value 
for spontaneous bending. This is what we are investigating in detail in this 
paper which is organised as follows. 

In section 2 we review the Mingaleev et al. [3] model . In section 3 
we study loop configurations in which the last two nodes of the chains are 
held together by a polaron. We extend this analysis to study hairpin-loop 
configurations in 4 for which the two opposite ends of the chain run parallel 
to each other, while the loop links the parallel strands together. Finally 
we show that one can estimate analytically the value of the parameters for 
which loops can be formed in section5. In the last two sections 6 and 7 we 
look in some detail at loop and hairpin-loop configurations for both single 
stranded DNA and polyacetylene and we show that the polarons, in these 
two systems, are very stable and that they can facilitate the formation of 
hairpin-loops. 
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2 Model 



The model proposed by Mingaleev [3] is described by the Hamiltonian 



U 1 " >R " ' + Un(R) ~ "Ul^l 4 +W I 2|0 n | 2 - J2 J, 



dr " v ' 2 



where i? n describes the position of each chain node, M is the node mass, W 
is the linear excitation transfer energy and A the non-linear self-trapping 
interaction. The excitation transfer coefficients J n ^ m are of the form: 

Jn, m = J{\Rn ~R m \) = (e a - 1) e -«\^-Z™\/\ (2) 

where a~ l sets the relative length scale over which the interaction decreases, 
in units of d, where a is the rest distance between two adjacent sites. The 
function J„ im describes the long range interaction between the electron field 
at different lattice sites n and m; its value decreases exponentially with the 
distance between them. Notice that when a is large and \R n — R m \ a, 
this corresponds to a nearest neighbour interaction with J n . m ~ <5n, m ±i(l + 

a(l - l ^"~^ m| )). 

In our formulation of the model, the normalisation of the electron field 
is preserved i.e. 

Ei^i 2 = 1 - ( 3 ) 

n 

The phonon potential U n consists of three terms: 
U n (R) = ^-(\R n -R n ^\-a) 2 + k °" 



2 U n_11 ; 2[i-(e n /e max ) 2 ] 

6 
2 



+^ ^(d- |it-i?m|) 2 e(d- |i?„-i? m |), (4) 



where the Heaviside function is defined as Q(x) = 1 for x > 1 and 6(x) = 
for x < 1. 

The first two terms in U n describe the elastic and the bending energy of 
the chain respectively, a is the equilibrium separation between nodes and 9 n 
is the angle between R n — i? n „x and R n -\.\ — Rn- Finally 9 max 

is the largest 

angle allowed between adjacent links. 



(1) 
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The term proportional to 8 in U n , models hard-core repulsion between 
the atoms of the chain. 5 should always be larger than a and d will corre- 
spond to the minimum distance allowed between nodes. 

For convenience, the symbols denoted by an overhead carat sign e.g. 
M, a etc. correspond to physical variables carrying units and dimensions 
while those without it correspond to dimensionless variables and parameters 
described below, except H, A and W which are dimensional quantities. We 
also use the symbol R for position of the nodes in physical units and r in 
dimensionless units. First we define the time scale r = hA/W 2 and use the 
lattice spacing a as the length scale. We can then define the dimensionless 
time t, position r and coupling constant g as 



t = T - 



A 

9= W 



- 

a 



(5) 



In terms of these variables the Hamiltonian takes the form 

2 



H 



a 4^ 



M fdf n 



+g ( 2|</> n | 2 - ]T Jnm<f>* n <t>m ] " ^"|<A. 



9 \m ,4 
n| 



(6) 



where 



v.<?> = - 2 (ft -?„-!-«) + 2 p _ 

+^ X ( d ~ l ffl ~ r m |) 2 9((i - \r n - f m \) 



(7) 



with 



M = M 
k = k 



a 2 W 2 
h 2 A 
A 

1^2 



„a 2 A 



,a 2 A 
V 2 



(8) 



Writing r n = (xi, n , ^2,n, ^3,n, ) we can derive the equation of motion for 
Xi tn from the Hamiltonian (6): 



M 



dt 2 
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k m<k aXi ' n 

-2<t>n+ ^2 Jnm<f>m + g\<f>n\ 2 <f>n = 0, (9) 

where the force F(t) and the friction term T dxi jn /dt, were added by hand to 
incorporate thermal fluctuations and F(t) was chosen as a delta correlated 
white noise satisfying 

< F(0)F(s) >= 2Tk B T5(s) (10) 

where 

k B T = k B Tj- = k B TWg. (11) 

As the equation for Xj is expressed in units of W 2 /(Ad), we have T = 
To 2 /h. The friction coefficient T can be evaluated from T ~ 67r^i?o where 
fj, = O.OOlPa s for water, and up to 4 times that value for the cytoplasm, 
where Rq is the average radius of a single molecule of the lattice. Notice 
also that the electron field 4> n is coupled to the phonon field x^ n through 
the function J nm . 

In what follows we are primarily interested in stationary configurations. 
To compute such solutions numerically we choose an initial lattice configu- 
ration with a loop structure and localised the electron so that it overlapped 
with both tails of the loop. We achieved this by using an approximation 
for the polaron electron field and distributing it over a few nodes spread 
between the two ends of chain. This way the polaron was able to bind the 
loop extremities together. We then relax the electron field keeping the lat- 
tice configuration unchanged and then evolved the entire system with an 
absorption term until it relaxed to a static configuration. This was achieved 
by solving equation (9) without thermal noise. 

In all our simulations we started from a very small value of k, typically 
k = 0.005, so that the lattice offered very little resistance to bending. We 
then increased the value of k in small increments using the relaxed confor- 
mation obtained for the previous k value as the initial configuration. We 
then equilibrated the system for the new value of k. By repeating the pro- 
cedure for each value of g we have determined the critical value k cr it(g) up 
to which the given configuration can be sustained by the polaron. 

Unless otherwise stated, we have used the following parameter values: 
5 = 10000, u = 1000 and M = 0.5. For stationary solutions the mass term 
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does not affect the results and 6 was chosen so that the repulsion potential 
is close to that of a hard shell. Finally for all the computed configura- 
tions, nothing prevents the nodes from being very close to their equilibrium 
distance and hence we have selected a relatively large value for a to ap- 
proximate stiff cross-node links. Following Mingaleev et al, we have also 
considered mainly the case a = 2 and d = 0.6. Finally we have also consid- 
ered the effect of varying the values of these two parameters. 

To solve equation (9) we used a 4th order Runge-Kutta method with 
a time step dt = 0.0001 in dimensionless units. To compute static con- 
figurations, we took T = 0, i.e no thermal noise, setting T = 1 and then 
integrating equation (9) until the system relaxed to a stationary solution. 

To study the thermal stability of the configurations for DNA and poly- 
acetylene, we solved equation (9), taking T = 300-fT and estimated T from 
the radius of the molecules as described above. For those simulations we 
started from the static configuration for which we wanted to evaluate the 
stability and let the system thermalise itself. The time needed for this ther- 
malisation was always orders of magnitude smaller than the average life 
time of the configurations we considered and so we did not need to resort 
to sophisticated thermalisation procedure as we did in [7] . 

3 Plain Loop Configurations 

Our first investigation involved considering a simple loop configurations for 
which all the nodes lie more or less on a circle with the two end points close 
to each other (separated by a distance d). When k is very small, the favoured 
configuration is one similar to the one presented in figure l.a. In this figures, 
the electron probability density is represented by the colour of the node. A 
dark colour corresponds to a null value while a light value corresponds to 
a higher probability density. The node at which the polaron field has its 
maximum value is close to the last 2 points at the opposite ends of the 
chain. This allows the electron field to be distributed on 2 nearby nodes 
rather than a single one and, as k is small, the deformation of the chain 
does not prevent this from taking place. As k increases, such a localisation 
becomes energetically expensive and the configuration assumes the shape of 
a horse-shoe, as presented in figures l.b. and I.e. As one increases k further, 
there is a point at which the stretching energy is too large and the polaron 
is not able to sustain the loop anymore. 

The difference between figures l.b. and l.c is that in the former the 
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(a) (b) (c) 

Figure 1: Loop configuration for N = 9 nodes for a) g = 2.5, k = 0.5, 
\M 2 = 0-407, \fo\ 2 = 0.221, |0 7 | 2 = 0.249; b) g = 2.5 k = 4, |0 O | 2 = |</> 8 | 2 = 
0.397; c) g = 5 k = 5, |^ | 2 = 0.642, \<p 8 \ 2 = 0.272. 

electron is localised equally on the 2 end nodes while in the later it is localised 
mostly on a single node. The difference is dictated by the value of g: for 
small <7, the polaron is wide and the electron spreads itself nearly equally 
between the two end points of the chain (Fig l.b). As g increases, the 
polaron becomes more localised and the electron becomes localised, more 
asymmetrically, on a single node (Fig l.c). 

The critical value of A; as a function of g is presented in Figure 2. a for 
loops consisting of 9 to 14 nodes. It is interesting to note that when g is 
small, the critical value of k is small. This can be explained by the fact that 
the coupling parameter g is small but also by the fact that the polaron is 
delocalised and hence the fraction of the electron close to the end point is 
smaller than for larger values of g. The maximum value of k cr u is reached 
for g ~ 5. For very large values of g, the electron is nearly fully localised on 
a single lattice point, but the attraction exerted by the polaron, surprisingly 
decreases very slowly. 

Having followed [3] and taken the values a = 2 and d = 0.6 for the re- 
sults presented so far, it is worth checking how these two parameters affect 
the results that we have obtained. We started by varying a, which controls, 
through J n>m , how fast the coupling between nodes decreases with the dis- 
tance separating them. The results are presented in Figure 2.b where we 
see that k cr u, contrary to what one might expect, increases with a. This is 
easily explained: having chosen d = 0.6, increasing a, not only reduces the 
long distance interaction between nodes but it also increases exponentially 
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the binding energy of nodes that are very close to each other. The binding 
energy of the end nodes, which are separated by a distance d < 1, thus 
increases with a. For this reason, we have decided to take d = a = 1 when 
we consider single stranded DNA and polyacetylene later in the paper. 

In figure 2.c we show how the critical value of k varies with d. As the 
parameter d sets the minimal distance allowed between 2 nodes and given 
that J n ^m decreases with the distance between nodes, it is not surprising 
that k cr n decreases when d becomes larger, but loop configurations can still 
be held by the polaron. 




g g g 



(a) (b) (c) 



Figure 2: Critical value of k for the existence of a loop configuration, a) 
a = 2 and N = 9 to U nodes, b) N = 9 nodes a = 1,1.5,2,25,3 and 
d = 0.6. 

c) N = 11 nodes d = 0.6, 0.75, 1 and a = 2. 



4 Hairpin-Loop Configurations 

Now we consider a hairpin-loop configuration as presented in Figure 3 similar 
to the structure that single stranded DNA can form and which is potentially 
more relevant to long chains. As for the plain loops, we generated these 
configurations for a small k and then slowly increased its value until the 
number of links, L, making the loop increased by one unit. This gave us the 
critical value, k cr n, for which the hairpin- loop configuration of a given size 
can be sustained by the polaron. 

The results are presented in figure 4 where we can see a sharp transition 
around g = 10. Below this value the hairpin-loop is only viable for relatively 
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small values of k but above it, they are sustainable for much more rigid 
chains. This is due to the fact that for small values of g, the polaron is 
always distributed over the handle of the hairpin-loop while when g > 10, it 
is localised mostly on one lattice site, at the base of the loop. In that case 
the interaction is stronger and supports loops for larger values of k. 

To make sure this was not an artefact of our procedure, we have tried to 
construct solutions using various initial conditions. We also used solutions 
obtained for g > 10 as initial conditions and then slowly decreased the value 
of g. Regardless of the procedure we used, we always obtained the curve of 
Figure 3. a. 

As expected, the configurations of figure 3 are harder to sustain than a 
simple loop as the chain needs to be bent near the stem of the hairpin-loop. 



Figure 3: Hairpin-loop configuration for N = 18 nodes. The brightness of 
the nodes is proportional to \<ft\ 2 . a) g = 1.5 and k = 0.5 (max \(p\ 2 = 0.188), 
b) g = 11 and k = 10 (max \<f>\ 2 = 0.887 ) 



5 Analytic Approximation 

Having computed numerically the critical value k cr n(g) for which the polaron 
is able to sustain a loop of a given size, we now try to estimate this value 
analytically. To do this, we consider a circular configuration of radius R 
made out of N segments, one of length b and N — 1 of length a, as depicted 
in figure 5. Note that the two nodes separated by the distance b are not 
linked to each other. If £ and [i are the angles opposite a and b, respectively, 





(b) 
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Figure 4: Critical value of k for the existence of an hairpin-loop config- 
uration, a) a = 2 and N = 7 to 11 node loops, b) N = 9 nodes and 
a = 1,1.5,2,2.5,3. 



we have 



(n-l)£ + /i = 2vr sin(|) = ^ sin(f ) = & (12) 



and so 



V = 6 ^ = 6 Sm ^ 2 ^ > 

Choosing specific values a = 1, b = 6q, this transcendental equation can 
be solved numerically to obtain the corresponding value £ = £o- We can then 
perform a first order expansions around this solution: b = bo + 5b, £ = £o + &; 
and obtain 

^ = — ? H & (W) 

^icos((iV-l)|)-|cos(|) 

To determine the critical value of g and k for which a loop configuration 
can exist, we have to minimise the Hamiltonian and, for each value of g and 
b, determine the value of k for which this Hamiltonian has a minimum. For 
each g, we then select the value of b for which k is the largest. 

Let us assume that the loop is symmetric, so that the N — 2 bending 
terms are all identical and are functions of £. The elastic terms are then 
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b 



Figure 5: Schematic representation of a polymer-loop with N nodes. N — 1 
bonds with rest distance a subtend an angle £ at the centre. The two end 
nodes, on which the electron is localised, are separated by a distance b and 
span an angle fi at the centre. 

also equal, but as they do not depend on £, they are constant and can thus 
be ignored for the minimisation. The repulsion term proportional to S can 
also be ignored if b > d. When b < d, the repulsion term leads to a very 
large energy increase and we can thus consider than b = d is the smallest 
value we should consider. 

To evaluate the electron field, we take the continuum limit of equation 
9 for stationary solutions: 



which is the well known non-linear Schoedinger equation, where g = g/{l — 
e~ a ), and which admits the following solution: 



Note that / \cf)(x)\ 2 dx = 1 and A = —g 2 /16. From the numerical solutions, 
we know that the wave function is centred on one of the two end lattice 
points and we can thus take (fro = </> c (0) and 4>n-i = <frc(bo)- 

The Hamiltonian can then be approximated by the following function of 



d 2 4>, 



■c 



+ g\(fic\ 2 4>c - \<Pc = 0, 



(15) 



dx 2 




(16) 
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[ht] N 


9 


10 


11 


12 


13 


14 


£o 


0.731341 


0.654972 


0.593071 


0.541875 


0.498824 


0.462115 



Table 1: Value of £o f° r various number of nodes N. 



b and £ 



H 



-2g(e a -l)e~ ab °/ a ^N-i + 



k (N-2)e 



1 



\ Umax J 



Next we compute the variation of H with respect to b and £ 



5H = 2ga</)o(t)N-i(e a - l)e- a0 °5b + 



\ u max 

Using equation (14) and imposing 5H = we get the condition 



(1-^(1 



1 



9 2 

"max 



(17) 



))«5£(18) 



k 



crit 



Jj<t>0<t>N-l 



(19) 



where 



D 



(N - 2)£o 



2a(e a - l)e~ ab ° 



\ Umax / 



(1-^(1 



e 2 

"max 



)) 



-4sin(|) 



(N — 1) cos((iV — 1)%) 



(20) 



2 , 6 cos(^) 

which depends on g but not on k. From equations (19) and (20) it is clear 
increases as 6o decreases and so we have to choose the smallest 
possible value for bo i.e. bo = d. 

Taking a = 1, bo = d = 0.6, we can solve equation (13) to obtain the 
values £o listed in Table 1. which we can use to estimate k cr u. The results 
are presented in figure 6 from which we see that our evaluation reproduces 
the gross features of the results obtained numerically (figure 4) : k cr ; lt is small 
when g is very small, and increases with g until a maximum is reached. Then, 
as g increases further k cr u slowly decreases. The maximum value obtained 
for k cr it is slightly smaller than the numerical value obtained before. The 
biggest discrepancy between the numerical and analytic results are for large 
g, but this is to be expected as this is the limit where the polaron is strongly 
localised and so is less well approximated by equation (16). 
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(a) (b) 

Figure 6: Theoretical estimation of the critical value of k for the existence 
of a loop configuration as a function of a and N. a) a = 2 and N = 9 to 14 
nodes, b) N = 9 nodes and a = 1, 1.5, 2, 25, 3 

6 Single Stranded DNA 

Having considered the Mingaleev et al. polaron model in general, we now 
consider two explicit cases: DNA and polyacetylene, both of which have 
parameters allowing the polaron to sustain loops. 

The parameters of our model were obtained from several sources. First 
of all, k can be determined from the flexural rigidity, k = AksT/Ro where 
Rq is the radius of the DNA strand. We do not have experimental values of 
5, but its actual value does not play an important role except that it must 
be large enough to mimic a hard shell repulsion. In practice, we chose a 
value larger than a. 

For single stranded DNA we have Rq ~ 0.33nm [8] A ss OAeV, W « 
0.3ey [3], k w O.lleV [8] and a « l.heV/A 2 [9]. This leads to the following 
dimensionless values: 

g « 1.33 a w 72.51 pa 0.487 M m 2.5 x 10 5 

k B T = k B T-£- « 0.115 T«3210 t « 2.92 x 10~ 15 s 

We thus see that single stranded DNA sits at the bottom left region 
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of figures 2, 4 and 6. For our simulations, we have chosen a = 2, 5 = 
100 and d = 1, the later parameter was taken as the worst case we could 
consider. We then found that DNA can easily sustain loops of 10 segments 
and hairpin-loops with 11 segments. We then studied the thermal stability of 
the configurations that we have obtained at T = 300-ftT. To achieve this, we 
started from a static configuration that we had obtained for DNA. We then 
solved equation (9), including the thermal noise, and ran 100 simulations 
for an extended period of time. We started by running 100 simulations for 
a loop made out of 10 nodes and we measured an average unfolding time 
of Ins. We also observed that a chain made out of 11 nodes is much more 
stable, and experiences a very slow unfolding of the loops with an average 
decay time of approximately 1.3/US. 

We have also performed thermal simulations for an hairpin-loop config- 
uration of 18 nodes and L = 11 and did not observe a single unfolding of 
the chain in 20/is. As the integration time steps was approximately 3 10~ 17 s 
this required 100 simultaneous simulations each performing around 10 10 in- 
tegration steps and we decided not to run it any longer as the stability of 
the loop was sufficiently well established. 

Under thermal noise, the stem, made out of the the two parallel ends 
of the chain, deforms itself and the chain takes the shape of a loop where 
the polaron links the two opposite ends of the chain around a couple of 
nodes, as presented in figure 7. In our simulations we observed that as the 
polaron moves along the chain, the size of the loop that it formed fluctuated 
constantly in time but it never unfolded. We can thus conclude that the 
DNA polarons loops are very stable. 

This loop configuration could play an important role in the formation of 
single stranded DNA hairpin-loops in vivo. The formation of such configu- 
rations depends on the likelihood of complementary DNA bases to face each 
other before they can bind by hydrogen bonding and this likelihood decreases 
rapidly as the length of the chain increases [10] and partially matching DNA 
base sequences are less stable than perfectly matched ones [11]. While homo- 
geneous sequences of DNA bases can bind quite rapidly like the one used in 
[10] and [11] for example, irregular sequences, like ATGCAGTC. . . GACTGCAT 
are less likely to match purely randomly. As the polaron folds the ssDNA 
into a loop and moves along the chain, the DNA bases on the opposite seg- 
ments of the loop slide relative to each other, under the action of thermal 
excitation, increasing the probability for a complementary sequence of bases 
to face each other and bind. In vitro, the reaction time of hairpin-loops has 
been determined to be several us [10], a length of time that, as we have 
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Figure 7: Thermalised, T = 300-fT, hairpin-loop DNA configuration for 
N = 18 nodes. The size of the disk are an exaggerated indication of their 
depth in the direction transverse to the plane of view. (The two nodes close 
to the crossing point are separated by the same distance as two neighbour 
points). The brightness of the nodes is proportional to \(f>\ 2 . 

shown, DNA-polaron can outlive easily. Hence, we can conclude that DNA 
polarons can increase the rate of formation of hairpin-loops. 

In our model, we have not taken into account the effect of water on 
the polaron. Recent studies [12, 13] have suggested that its effect would 
be to reduce the polaron size, which, as can be seen from equation (16), 
corresponds to increasing the value of g. The net effect would thus be to 
move DNA to a parameter region where the polaron is stronger, as can be 
seen in figures 2 and 4. As a result the polaron would then be able to sustain 
smaller loops. 

7 Polyacetylene 

For polyacetylene, the physical value of the parameters are given by Ro ~ 
0.24nm, W « 2.5eV, a rj 2leV/A 2 [14], k » 3.7eV/ [15] and so 

g ss 2.56 a « 128 k ~ 3.79 M w 3.2 x 10 4 

k B T w 0.026 T w 1316 r f» 6.74 x 10~ 16 s 
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Once again, we took d = 1 and a = 2 to avoid the potentially spurious 
effects induced at close distance by J n ,m- In this case, we were able to make 
loops out of 12 nodes and hairpin-loops of 12 nodes too. 

Under thermal fluctuation, both were very stable. In this case, the inte- 
gration time step was approximately 710~ 18 s and our attempt to evaluate 
the configuration average life time was achieved by running 100 simultaneous 
each performing over 10 10 integration steps and we did not observe a single 
unfolding of the chain in 13fis. We also ran simulations for an hairpin-loop 
configuration of N = 18 nodes and L = 12 and also did not observe a single 
unfolding of the chain in over 10/US. 

8 Conclusions 

In this paper we have studied the possibility of a polaron to sustain loops 
and hairpin-loop configurations. In these configurations the polaron was 
localised over lattice nodes that were well separated along the chain back- 
bone but spatially close to each other because of the bending of the chain. 
The polaron then acted as a linker between the two regions of the chain 
and so could sustain the loop configuration if the chain was not too rigid. 
The Mingaleev model we have used to describe this property takes into ac- 
count the long distance interactions between the electron and the phonon 
field, with a strength decreasing with the distance. For the configurations 
we have studied, the most important contribution comes from lattice nodes 
that are spatially close to each other and the energy contribution from next 
to nearest neighbour is not essential, unlike in our study of spontaneous 
polaron displacements [7] where the next to nearest neighbour terms were 
essential for the polaron to move along the bending gradient of the chain. 

We have determined the critical value of the chain rigidity k cr n as a 
function of the polaron coupling constant g and we have shown that polarons 
are able to sustain relatively small loops for a wide range of parameters 
values. 

We have then shown that both DNA and polyacetylene are flexible 
enough for a polaron to sustain hairpin-loop configurations. Moreover, we 
have also shown that these hairpin-loop configurations are very stable under 
thermal excitations, with average live times exceeding lO^us, and that they 
can facilitate the formation of hairpin-loops of single stranded DNA. 
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